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Abstract 


A detailed theoretical study of the optical absorption in self-assembled 
quantum dots is presented in this paper. A rigorous atomistic strain 
model as well as a sophisticated electronic band structure model are 
used to ensure accurate prediction of the optical transitions in these 
devices . The optimized models presented in this paper are able to 
reproduce the experimental results with an error less than 1%. The 
effects of incident light polarization, alloy mole fraction, quantum dot 
dimensions, and doping have been investigated. The in-plane polar¬ 
ized light absorption is more significant than the perpendicularly po¬ 
larized light absorption. Increasing the mole fraction of the strain 
controlling layer leads to a lower energy gap and larger absorption 
wavelength. Surprisingly, the absorption wavelength is highly sensi¬ 
tive to changes in the dot diameter, but almost insensitive to changes 
in the dot height. This unpredicted behavior is explained by sensi¬ 
tivity analysis of different factors which affect the optical transition 
energy. 

Keywords — Self-assembled quantum dots, Stranski-Krastanov, Anhar- 
monic atomistic strain model, Biaxial strain ratio, Semi-empirical tight bind¬ 
ing, sp3d5s* with spin-orbit coupling (sp3d5s*_SO), Quantum dot hlling. 
Optical absorption. 
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1 Introduction 


Self-assembled quantum dots have improved the performance of many op¬ 
toelectronic devices, such as quantum dot infrared photodetectors (QDIPs) 
(^[^, and intermediate band solar cells (IBSCs) [^|^. QDIPs have lower dark 
current than conventional photodetectors and are sensitive to normally 
incident light unlike their counterparts that are made from quantum wells . 
For IBSCs, one of the most successful methods in pushing solar cell efficiency 
beyond the Shockley-Queisser limit is to add one or more intermediate bands 
inside the gap, which can be realized by using quantum dots and quantum 
dot-in-a-well devices [^. 

The absorption coefficient a(A) of quantum dots is an important param¬ 
eter that needs to be well designed for the proper operation of these devices, 
and an accurate model for the absorption coefficient a(A) would be highly 
appreciated by researchers as well as engineers working in these helds. For 
this reason, this study has been devoted to £11 in some of the gaps found 
in the available models of the absorption, specially in the atomistic strain 
calculations that are needed for accurate description of the bound states. 
In addition, a detailed study is provided of the effects of varying doping, 
dimensions, alloy mole fraction, and incident light polarization. 

Self-assembled quantum dots are highly strained heterostructures, the 
atomistic strain in such structures is usually on the order of 10%. Such high 
strain values are beyond the domain of validity of the continuum elasticity 
theory . Effects of interface inter-diffusion and alloy disorder further com- 
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plicate the physics of determining the relaxed atom positions and a rigorous 
atomistic treatment of strain is needed (such as Keating or anharmonic 
models 0 )- 

The anharmonic strain model reported by Lazarenkova, et. ah mod¬ 
ified the standard harmonic Keating model to include the effect of anhar- 
monicity in the lattice potential. This modification is necessary as the har¬ 
monic Keating potential underestimates the repulsive forces at smaller atomic 
separations and fails to reproduce the weakening of the strain energy at large 
bond lengths [^. The anharmonic strain model of Lazarenkova et. ah has 
hve strain parameters for optimization compared to just two parameters in 
the harmonic Keating model. The anharmonic strain parameters were orig¬ 
inally optimized to obtain correct Griineisen parameters for more accurate 
phonon dispersion calculations [^, however, it will be shown later that the 
original parameter set cannot reproduce the experimental optical absorption 
peaks in quantum dots [^. Moreover, these parameters do not capture the 
correct strain distribution in quantum wells, for which an analytical solution 
is known. The approach described in this paper for optimizing the parameters 
to obtain correct biaxial strain ratio in the quantum wells has resulted in an 
improvement in quantum dot simulations in terms of capturing experimental 
optical transitions much more accurately than previously. The Hamilto¬ 
nian has been constructed with semi-empirical tight binding with 20-orbital 
sp3d5s* basis per atom, including spin-orbit interaction (sp3d5s*_SO) [lo] . 
The absorption coefficient has been calculated by employing Fermi’s golden 
rule. 


4 


In section II, the systems simulated and the numerical tools used are 
described. In section III, the theoretical aspects of the problem and the 
optimization of the strain model is discussed. Lastly in section IV, the results 
of the simulation is discussed in addition to analyzing the sensitivity of the 
absorption to various quantum dot parameters. 


2 Multi Million Atom Simulation 


As shown in the Figure the investigated SK-QD system 11 has a dome 
shaped InAs quantum dot with a base diameter of 20 nm and a height of 
5 nm. The wetting layer is 2 monolayers. The measured system has been 
doped with sheet doping of two electrons per dot. The strain controlling layer 
is made of Ino.isGao.ssAs and is sandwiched between two layers of GaAs each 
with a thickness 1 nm. Next, there are two layers of Alo. 22 Gao. 78 As, each with 
a thickness of 2 nm. The rest of the structure is made of Alo. 07 Gao. 93 As. The 
dimensions of the simulated SK-QD systems are 60 nm x 60 nm x 60 nm. 
The strain simulation contains around 10 Million atoms and the atomistic 
grid is shown in Figure]^ The band structure calculations do not need all of 
these atoms, since bound states decay exponentially outside the quantum dot. 
The band structure calculations are performed on 40 nm x 40 nm x 20 nm 
box surrounding the quantum dot, this box contains only 1.5 million atoms. 
Each atom has 20 orbitals in the sp3d5s*_SO tight binding basis. Such large 
systems are computationally expensive to run strain and electronic structure 
simulations and require highly scalable computational codes, the described 
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system would take around three hours on 729 cores. The code that has 
been used for our simulations is the Nano Electronic Modeling tool version 


5 "NEMOS” 
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3 Theoretical Model 


3.1 Atomistic strain model 


The Harmonic Keating strain model treats the interatomic forces as spring 
forces connecting the atoms together. It considers nearest neighbor interac¬ 
tions only and the expression for the elastic energy is given by 




^mn ( 2 
"J2 Vmn 

in 


~ ^mn) + X] 


k>n 


f^mnk 

dmndmk 


(^mn ‘ ‘^mk ' ^m/c) 


, ( 1 ) 


where the coefficient a corresponds to the force constant of the bond length 
distortion, and {3 corresponds to the bond angle distortion as shown in Figure 
1 ^ a and (3 are material-dependent constants, r^n is the displacement vector 
from atom m to atom n for the strained crystal, while dmn is the same 
vector for the bulk unstrained crystal. The summation is taken over the 
nearest neighbors only and the total energy is minimized with respect to 
the individual atomic positions, relaxing the structure. The problem with 
the harmonic Keating potential given by equation ([^ is that, it produces 
a symmetric energy prohle around the equilibrium interatomic distance and 
angle. For this reason, the Keating model fails to reproduce the weakening 
of the strain energy with increasing bond length and it underestimates the 
repulsive forces at close atomic separation The anharmonic correction 
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of the Keating model proposed by Lazarenkova, et. al. has been able to 
solve this problem by modifying the two parameters a and (3 of the Keating 
model and making them functions of bond length r and bond angle 9 as given 
by 

/ 

a = \ 1- A — _^ 

^mn ^mn \ ^ -^mn 12 

\ ^mn 

A™* = (1 - (cos («„„0 - cos (CJ)) X (1 - 

\ ^mn ^mk 

In this anharmonic strain model, there are hve strain parameters for each 
material, B^nk, and Cmnk- The anharmonic model was 

developed mainly for simulating phonon dispersion and transport. For this 
reason, the anharmonic strain parameters were optimized to reproduce the 
Griineisen parameters 7 * = which are a measure of the dependence 

of the phonon mode frequencies on strain. Simulating the strain in quantum 
dots with the original anharmonic strain parameters provide inaccurate re¬ 
sults that do not match experimental results as shown in the results section 
of the paper. 

In addition, simulating strain in quantum wells with these parameters 
gives strain tensor components that do not match well the analytical solution 
of the strain in quantum wells as shown in Table The parameters of the 
model described in this paper have been tuned to reproduce the biaxial strain 
ratio a of InAs in order to capture the strain distribution in quantum wells 
and quantum dots. The InAs biaxial strain ratio a of 1.053 was obtained 
from more accurate but size-limited density functional theory calculations 
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Method 

^11 

e± 

Analytical 

- 6.68 % 

7.04 % 

Anharmonic before tuning 

- 6.68 % 

8.9 % 

Anharmonic after tuning 

- 6.68 % 

7.04 % 


Table 1: Strain calculated for InAs/GaAs quantum well. Tuning has im¬ 
proved the anharmonic strain results in quantum well. 


13 . Only one parameter needs to be optimized in this approach. Out of 


the hve strain parameters, tuning the parameter while keeping the other 
parameters as reported in yields the best convergence and stability of the 
strain model. The new value of after optimization is 19.35Nm“^, which can 
be readily used in accurate large-scale atomistic strain simulations. Table 
shows the atomistic strain calculated for InAs/GaAs quantum well, as shown 
in Figure before and after tuning. The analytical expressions for the strain 
components in quantum wells are en = °‘c:aAs-a.i„As _ —(j^n 14 ^ where 

^InAs 

a is the lattice constant. 


3.2 Electronic structure and absorption 


The eigenstates of the system were calculated with a Hamiltonian constructed 
with semi-empirical tight binding sp3d5s*_SO basis. The Slater-Koster TB 


15 parameters for InAs, GaAs, and AlAs are taken from 16,17. For includ¬ 


ing the effect of strain on the tight binding Hamiltonian, please refer to 10 


These parameters are well established and have been verihed before versus 
















experimental measurements of quantum dots 18 -21 


For the absorption coefficient a, Fermi golden rule has been used to cal¬ 


culate the absorption coefficient 22,23 
2 TT Ca 


a(uj} = 


n eo c 


T.\i 


fi ■ h 


S{Ef-E,-hio){Fi-Ff) , 


( 3 ) 


where ridots is the number of quantum dots per unit volume, oj is the 
photon angular frequency, Ei and Ef are initial and final energies of the 
transition, Fi and Ff are occupation probability of the initial and final states, 
n is the refractive index of the material, c is the speed of light in free space, 
Co is the free space permittivity, e is the polarization of the incident light, 
and dfi is the first order dipole moment that is given by d/* = q < >, 

where q is the electron charge. 

For the transitions between the valence bound states and the conduction 
bound states, Fj = 1, where Ff depends on the energy level and the doping. 
Normally, quantum dots are occupied by the number of electrons equal to 


the average number of dopants per dot 24 . This approach is reasonable 
for quantum dots that are far from heavily doped regions, however, it is 
not appropriate for quantum dots adjacent to heavily doped regions, such 
as contacts. In addition, it is assumed that if one level is hlled then the 
next electron fills a higher level instead of filling the same level with opposite 
spin due to the high Coulomb repulsion between them. Irrespective of the 
Filing, the Coulomb repulsion shift in the energy has been neglected. In other 
words if the doping is two electrons per dot, then the electron ground state 
and the first excited state will be occupied each by an electron and they will 
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not accept more electrons. These assumptions have been considered while 
calculating the absorption. 

4 Results and Discussion 

4.1 Simulation versus Experimental Results 

A comparison with the measured absorption spectrum of the SK-QD 
system has been made to validate the model. Figure shows the calculated 
and measured absorption spectrum of the device. The simulation result 
matches very well with the measured absorption and the error in estimating 
the energy of the absorption peak is less than 1 %. Further comparison with 
experimental measurements are provided later when discussing the effect of 
the alloy mole fraction of the strain controlling layer. 

4.2 Band structure and states 

Figure shows the wavefunction probability density of the hrst eight non¬ 
degenerate states of both electrons and holes. It is worth noting that the 
hole ground state has an s orbital like shape. 

SK-QDs have complicated band profile as multiple effects, such as ge¬ 
ometric confinement, strain, alloy disorder, etc. can cause major changes 
in the band edges of the bulk material. It is important to know where the 
wavefunctions of the electrons and holes are localized due to these disordered 
band edges, as the spatial overlap between the states will determine the opti- 
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cal absorption spectrum. Hence, one can look at the conduction and valence 
band edges along arbitrary lines passing through the quantum dot. This can 
be done by using the deformation potential theory that gives the shift of the 
band edges due to small lattice deformations. The shift in the band edges 


due to lattice strain for zincblende materials is given by reference 25 

/\Ec ) 


^E^HH — Clv^H + 2^B, 


(4) 


^E^lh — o,vBh — 


where AEc is the shift in the conduction band edge, AE^hh and AEyin 
are the shifts in the heavy and light hole band edges, respectively. Oc, a„, 
and b are the deformation potential coefficients of the material. In these 
simulations, we have used the parameters recommended for III-V materials 


by 26 . Eh and Eb are the hydrostatic and biaxial strain components which 
are linear combinations of the atomistic strain components: Eh = £xx + Byy + 
Ezz and Eb = £xx + Byy — 2 Ezz 25 , where z is the growth direction. Figure]^ 
shows the band edges along two lines through the middle of the quantum dot 
along the [001] and [110] directions. The unstrained band edges are plotted 
to show the signihcant effect of strain on the band edges. 


4.3 Doping and Polarization 

Figure [^shows the absorption spectrum of the SK-QD with different cases of 
doping for in-plane polarized light (normally incident light) and perpendicu¬ 
larly polarized light. The in-plane polarized light absorption is more than ten 
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times larger than the perpendicularly polarized light absorption. Quantum 
dots can absorb both polarization components, but their main advantage 
over quantum wells are the dots’ sensitivity to normally incident light which 
makes them more efficient than quantum wells in applications like solar cells 
and photodetectors. For this reason, we will focus more on studying the in¬ 
plane polarized absorption properties. In Figure]^ each peak is coming from 
transitions between valence states and each conduction state. Practically, 
these peaks can overlap with each others if the inhomogeneous broadening, 
due to variance in the dot dimensions, is high. 

4.4 Effect of Quantum Dot Dimensions 

Figure shows the effect of changing quantum dot diameter and height on 
the in-plane polarized absorption spectrum. Increasing the dot diameter 
results in red-shifting the peaks, while increasing the dot height doesn’t have 
a signihcant effect on wavelength. In contrast, to the simple particle in a box 
problem, which predicts a stronger sensitivity to the smaller dimension (the 
height), this work shows that absorption wavelength is much more sensitive 
to changing the dot diameter than the dot height. 

The effects of changing the dimensions on the energy transition AE be¬ 
tween the hole and the electron ground states can be understood with a 
simple analytical model. This transition has two contributions: strain and 
conhnement. The strain shifts the band edges and affects the energy gap Eg 
, while the conhnement increases the minimum allowed energy of electron 
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Eelec and hole Ehoie with respect to band edges. Let Ebox = E^iec + Ehoie, 
then the transition energy E is 

/S.E = Eg + Ebox- (5) 

Due to the sign of the deformation potential and strain, the valence band 
edge inside the quantum dot is the heavy hole , from equation (|^ 

b 

Eg (^Ctc 2 '^^* 

Figure [TO] shows the effect of changing the diameter and the height on the 
hydrostatic and biaxial strain. The magnitude of the biaxial strain increases 
with increasing the diameter and decreases with increasing the height, while 
the magnitude of the hydrostatic strain changes slightly in the opposite di¬ 
rection to the biaxial strain. Increasing the height is equivalent to decreasing 
the diameter in terms of changing the strain in the quantum dot. Increas¬ 
ing the diameter reduces the energy gap which further reduces the optical 
transition energy, while increasing the height increases the energy gap which 
works against the reduction in the conhnement energy. This compensation 
results in almost the same optical transition energy. 

Although the variations in the hydrostatic strain are smaller than the 
variations in the biaxial strain, as shown in Figure [TOl the hydrostatic strain 
variations shouldn’t be neglected. This is due to the higher deformation 
potential weight for the hydrostatic strain. For example, Oc — = —6eV 

is six times higher than | = —leV for InAs. Also, changing one of the 
dimensions either increases or decreases the hydrostatic strain, and it will 
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have the opposite effect on the biaxial strain (decreases or increases), but 
the hydrostatic and biaxial strain will work together in the same direction 
on the energy gap since their terms have opposite sign in equation ([^. 

To get an expression for the Ehox, the dome shaped quantum dot is ap¬ 
proximated to be a disc of cylinder R and height H, one can easily obtain 
Ebox by solving an effective mass Hamiltonian in the cylindrical coordinates, 


( 1 


Ebox — vr-1-+ 


rrip 


rrih 


TT 


X, 


01 




rrih 


(7) 


where rrie, rrih are the electron and the heavy hole effective masses, and 
Xoi = 2.405 is the hrst zero of Bessel function of the hrst kind with order 
0. It should be noted that effective masses for the electron and heavy hole 
under strain are different from the bulk, and the InAs effective masses are 


nie = O.lrrio and rrih = 0.48mo 16 


For a quantum dot of H = 20 nm and H = 5 nm, the sensitivity of the 


conhnement energy to the dot radius is 


SE. 


SR 




mh 




meV/nm, and the sensitivity of the energy gap to the dot radius is ^ = 
(oc—a.y)^ —~ —17 meV/nm which give a total sensitivity of the optical 
transition —22 meV/nm. As for changing the height, in the same way 

— = + -68 meV/nm, ^ = (a, - ~ 


5H 


me ' m^ J 

5AE 


65 meV/nm which give ~ —3 meV/nm. For increasing the dot radius, 
both contributions reduce the transition energy. While increasing the dot 
height, both contributions are working against each other which reduces the 
sensitivity of the transition energy to the dot height. 
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4.5 Strain Controlling Layer 


Changing the In mole fraction of the InGaAs strain controlling layer (capping 
layer) is a convenient way to tune the absorption peak. The effect of mole 
fraction has been studied on a slightly different system, reported in ref. [^, 
which help us further validate the results of the simulations. The system 


reported in ref. 9 is almost the same as the system in 11 except for two 


differences: First, it does not have any doping. Second, instead of AlGaAs 
alloys we have only GaAs material. Figure [TT] shows the experimental and 
simulation results of the optical transition of the SK-QD system reported 
in ref. [^. The optimization of the a n h armonic strain model has greatly 
improved the simulation results. Increasing the In mole fraction increases 
the transition wavelength. To understand the reason behind this, one needs 
to see the effect of changing the mole fraction on the hydrostatic and biaxial 
strain and hence the band edges. Figure shows the hydrostatic and biaxial 
strain along two lines passing through the middle of the quantum dot in the 
[001] and [110] directions for different cases of In mole fractions. As shown 
in these figures, the hydrostatic and biaxial strain change with the In mole 
fraction in the same way they change with diameter; increasing the In mole 
fraction results in an increase in the magnitude of the biaxial strain and a 
decrease in the magnitude of the hydrostatic strain. This leads to a lower 


energy gap and larger absorption wavelength, as shown in Figure 11 
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5 Conclusion 


In this paper, a detailed theoretical study of the optical absorption and strain 
behavior in self-assembled quantum dots has been presented. Self-assembled 
quantum dots are highly strained heterostructures, and a rigorous atomistic 
strain model is needed to accurately calculate the electronic states in the sys¬ 
tem. We have described an optimization procedure of the anharmonic strain 
model, which has greatly improved the capability of simulations to reproduce 
experimental optical characteristics. The simulation is able to reproduce the 
experimental results with an error less than 1%. The optimized strain model 
has been implemented in NEMOS and has been used to simulate character¬ 
istics of an InAs/GaAs/AlAs quantum dot systems. The inplane-polarized 
light absorption is more signihcant than the perpendicularly polarized light 
absorption. Increasing the dot diameter results in shifting the peaks towards 
longer wavelengths, while increasing the dot height doesn’t seem to have a 
signihcant effect on wavelength. In case of changing the diameter, changes 
in band gap and conhnement energies work with each others, while in case 
of changing the height, the changes in band gap and conhnement energies 
work against each others. Increasing the In mole fraction in the strain con¬ 
trolling layer works in the same way as increasing the dot diameter in terms 
of changing the strain which leads to longer absorption wavelengths. 

In conclusion, the method presented here provides a way to incorporate 
the inhomogeneous environment of QDs in simulations by taking into account 
device geometry and quantum conhnement, alloy disorder, electrostatics, and 
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spatially varying strain distribution. Such details are needed to interpret 
and guide experimental measurements and device design with quantitative 
accuracy. 
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Figure 1: A schematic of the measured and simulated SK-QD system. The 
dimension of the simulated structure is 60 nm X 60 nm X 60 nm. The 
quantum dot is a dome shaped InAs with base diameter of 20 nm and height 
of 5 nm, while the wetting layer is 2 monolayers. The strain controlling 
layer of Ino. 15 Gao. 85 As is sandwiched between two 1 nm layers of GaAs, and 
two 2 nm layers of Alo. 22 Gao. 78 As. The rest of the structure is made of 
Alo. 07 Gao. 93 As. 


23 










n 



^mk 


m 


Figure 2: Bond lengths and bond angle for three neighboring atoms m, n, 
and k. 



Figure 3: An InAs/GaAs Quantum well of thickness 3 nm used for the 
optimization of the anharmonic strain model. . 
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Figure 4: The atomistic grid of the simulated SKQD showing the different 
random alloy regions. The number of atoms used for strain simulation 
10 Million atoms, while for band structure calculations only 1.5 Million 
atoms are needed. Lines LI and L2 are two lines passing through the middle 
of the quantum dot in the directions [001] and [110] respectively. 
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1.2 



Figure 5: This figure shows the calculated and measured 11 absorption of 
the SK-QD system. The quantum dot is dome shaped with base diameter of 
20 nm and height of 5 nm. The doping is 2 electrons per dot. The calculated 
absorption matches well with the experiment, the error is less than 1% in the 
absorption peak. 
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Figure 6: The conduction and valence band edges (solid lines) along a line 
through the middle of the quantum dot in the [001] and [110] directions. The 
dashed lines are the band edges of the unstrained bulk materials, drawn to 
show the significant effect of strain on deforming the band structure. Also 
shown are the lowest energy bound states in the quantum dot. 
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Figure 8: The absorption spectrum at different cases of doping for in-plane 
(A) and perpendicularly (B) polarized incident light. 
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Figure 9: The in-plane polarized absorption spectrum calculated for different 
diameters (A) and different heights (B) of the quantum dot. Increasing 
the dot diameter results in shifting the peaks towards longer wavelengths, 
while increasing the dot height doesn’t seem to have a significant effect on 
wavelength. 
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Figure 10: Hydrostatic Eh and Biaxial Eb strain with different dimensions 
along a line through the middle of the quantum dot in the [001] direction. 
Figure A is Eh at different dot diameters, Figure B is Eh at different dot 
heights, Figure C is Eb at different dot diameters, and Figure D is at 
different dot heights. The magnitude of the biaxial strain increases with 
increasing the diameter and decreases with increasing the height, while the 
hydrostatic strain behaves in the opposite way. 
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Figure 11: Experimental and simulation results of the optical transition of the 
SK-QD system reported in ref. [^. Increasing the In mole fraction increases 
the transition wavelength. The optimization of the anharmonic strain model 
has greatly improved the simulation results. 
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Figure 12: Hydrostatic and biaxial strain with different In mole fraction 
along two line through the middle of the quantum dot in the [001] and [110] 
direction. 


Figure captions 

• Figure A schematic of the measured and simulated SK-QD system. 
The dimension of the simulated structure is 60 nm X 60 nm X 60 nm. 
The quantum dot is a dome shaped InAs with base diameter of 20 
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nm and height of 5 nm, while the wetting layer is 2 monolayers. The 
strain controlling layer of Ino.isGao.ssAs is sandwiched between two 1 
nm layers of GaAs, and two 2 nm layers of Alo, 22 Gao. 78 As. The rest of 
the structure is made of Alo.orGao.gsAs. 

• Figure Bond lengths and bond angle for three neighboring atoms m, 
n, and k. 

• Figure]^ An InAs/GaAs Quantum well of thickness 3 nm used for the 
optimization of the anharmonic strain model. 

• Figure The atomistic grid of the simulated SKQD showing the dif¬ 
ferent random alloy regions. Lines LI and L2 are two lines passing 
through the middle of the quantum dot in the directions [001] and 
[110] respectively. 

• Figure This hgure shows the calculated and measured absorption 
of the SK-QD system. The quantum dot is dome shaped with base 
diameter of 20 nm and height of 5 nm. The doping is 2 electrons per 
dot. The calculated absorption matches well with the experiment, the 
error is less than 1% in the absorption peak. 

• Figure The conduction and valence band edges (solid lines) along 
a line through the middle of the quantum dot in the [001] and [110] 
directions. The dashed lines are the band edges of the unstrained bulk 
materials, drawn to show the signihcant effect of strain on deforming 
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the band structure. Also shown are the lowest energy bound states in 
the quantum dot. 

• Figure The magnitude square of the wave functions of the electron 
and hole state. Plotting only the hrst eight electron and hole states. 

• Figure The absorption spectrum at different cases of doping for 
in-plane (A) and perpendicularly (B) polarized incident light. 

• Figure The in-plane polarized absorption spectrum calculated for 
different diameters (A) and different heights (B) of the quantum dot. 
Increasing the dot diameter results in shifting the peaks towards longer 
wavelengths, while increasing the dot height doesn’t seem to have a 
signihcant effect on wavelength. 


Figure IT Hydrostatic Eh and Biaxial Eb strain with different dimen¬ 
sions along a line through the middle of the quantum dot in the [001] 
direction. Figure A is Eh at different dot diameters, Figure B is Eh at 
different dot heights. Figure C is at different dot diameters, and Fig¬ 
ure D is at different dot heights. The magnitude of the biaxial strain 
increases with increasing the diameter and decreases with increasing the 
height, while the hydrostatic strain behaves in the opposite way. 


• Figure Experimental and simulation results of the optical transi¬ 
tion of the SK-QD system reported in ref. [^. Increasing the In mole 
fraction increases the transition wavelength. The optimization of the 
anharmonic strain model has greatly improved the simulation results. 
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• Figure [T^ Hydrostatic and biaxial strain with different In mole fraction 
along two line through the middle of the quantum dot in the [001] and 
[110] direction. 
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